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dynamics of the system. We show that the long time tail of the decay is 

algebraic (1/i), irrespective of the dimension D. 
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I. INTRODUCTION 

In a previous work [|l| we studied the decay of quasibounded classical hamiltonian systems 
in two dimensions, in particular the decay problem for the Sinai well. In this report we extend 
the study to the D-dimensional case. We will briefly describe what a quasibounded system 
means, but for further details we refer to the original work. 

A quasibounded system is a dynamical system transiently bounded to a finite region of the 
phase space where an infinite set of non stable periodic orbits is included before it displays 
unbounded dynamics. The transition from the bounded motion to the unbounded one is the 
decay process, and the decay law is related to the bounded transient dynamics. The kind of 
systems that we are interested in is fully chaotic but not completely hyperbolic. In terms of 
the invariant set it means that instead of being completely hyperbolic, it has a non hyperbolic 
(namely parabolic ||) subset. One of the main differences with the analogous system in two 
dimensions is that whereas in that case the invariant set could be fully hyperbolic or have a 
parabolic subset depending on the value of a simple parameter, for the D > 2 dimensional 
system the invariant set has always a parabolic subset. In this case the global consequence 
on the decay law is that it always exhibits a crossover between a stretched exponential and 
an algebraic decay for long times. One of the main purposes of this work is to study the 
long time tail of the decay law in D dimensions and to relate the decay of population from 
equilibrium to internal distributions that characterize the dynamics. 

Our work is organized as follows. In Sec. |T| we introduce the system that we will call the 
Sinai well in D dimensions because it is geometrically similar to the Sinai billiard |§ in D 



dimensions but with a finite rather than infinite well. In Sec. [Ifl] we review some results of 
the analogous model in two dimensions fLf. 



Sec. |Iy| is the central body of the work and is devoted to the study of the temporal decay 
law both numerical and analytically. We show the results of the numerical study of the 
decay that reveal the above mentioned behavior, and using ergodic properties we relate the 
decay law to internal distributions that depend on the internal dynamics. Finally, Sec. [V] 



is devoted to discussions and conclusions. We include Appendix A in which the explicit 
expression for the transition probability from the bounded to the free region is derived. 

II. THE SINAI WELL IN D DIMENSIONS 

Our system consists of a point particle of unit mass in a D (D > 2) dimensional square 
well (hypercube) of depth — Vq, (Vq > 0) and side a, which collides elastically with a fixed 
scatterer located in the center of the well. The scatterer is a D dimensional unit sphere and 
as usual we take the speed of the particle to be one. In the study of the analogous system 
for D = 2 [p]] we consider the total energy E = p 2 /2 — Vq (0 < E < Vo), and explain in detail 
how the collision with the scatterer rearranges the energy among the degrees of freedom in 
order that the particle could transit from the bounded region to the free one after colliding 
with the central barrier. 

Although the extension to higher dimensions is straightforward, there are some new features 
to remark on and we devote the remainder of this section to these results. 
The bounded motion in D dimensions is caracterized by the condition 



1_ fP_^l\ >sin^ m , (2.1) 

\ V |p| / 

where rij, i = 1, .., D , is the normal direction to face i on which the particle bounces and 

^um = arcsin , (2.2) 

j l + V /E 



is the limit angle in D dimensions. The inequality (|2.1|) is the condition to have an internal 
reflection when the particle reaches the boundary of the hypercube. This result has been 
explained in detail for D = 2 in Appendix A of 0, and for D > 2 the derivation is similar. 
As we mentioned before, the collision with the scatterer changes the value of the components 
of the momentum and this enables the transition to the free region, or in other words the 
decay process. In the D dimensional problem the limit angle ^u m can be related to the 
probability ujd < 1 that the particle transits from the bounded to the free region after one 
collision with the scatterer. The space of momenta is a D dimensional unit sphere in which 



UD = ~im~ ' ( 3) 

where d^u m (D) and Q(D) are respectively the solid angle subtended by ^u m and the total 
solid angle in D dimensions. In Appendix A we derive an explicit expression for (|2.3|) in 
terms of ^u m . Here we give the result 

"d ~ ^L 1 (2-4) 

which shows that the transition probability decreases with the dimensions D when 
$(i m <C 1 for a given (fixed) energy E. 

III. PRELIMINARY REMARKS 

In this subsection we review some results of the analogous problem in two dimensions [|l| . 
The decay of population N(t) inside the well is characterized by two distributions specific 
to the internal dynamics. The first one g(t)dt is the fraction of particles for which the first 
collision with the circular scatterer occurs between t and t+dt, and the second one f(t)dt is 
the fraction of particles for which the time between two successive collisions with the central 
scatterer lies between t and t+dt. In the decay process from the equilibrium population 
these distributions are not independent, but related by 

dg/dt = -g(0)f(t) , (3.1) 

In UU (in the following I) we conclude that distributions g(t)dt that decrease exponen- 
tially or faster lead to an exponential decay law while an algebraic decay of g(t) gives rise 
to an exponential decay law that changes into a power law decay for long times. Our model 
is closely related to the periodic Lorentz gas in two dimensions, in which the problem of the 
asymptotic behavior of the velocity autocorrelation function has been studied extensively 
over the last several years both theoretically and numerically (see 0] and references cited 
there). In that model we can distinguish two kinds of behaviors depending on the "finite" 
or "infinite" nature of the horizon. By definition a periodic configuration of scatterers has 



infinite horizon if the length of the free motion of the particle is unbounded. Actually if 
the horizon is infinite then there exist trajectories that never reflect from the scatterers. 
These trajectories define the so called corridors that are characterized by the directions of 
the velocity such that v y /v x = Z\jz% (here z\ and z-i are coprime integer numbers). The 
number of open corridors increases when the radii of the scatterers R decrease. Following 
the notation of ref. [||] for the square Lorentz model there are at least two of such open 
channels corresponding to x-direction and y-direction, and we call these channels a and (3 
respectively. When a/5/10 < R < \/2/4 we have other open channels (7) that correspond to 
v y /v x = ±1. In I we explained in detail the connection between the open corridors and the 
parabolic periodic orbits that appear in our system when the radius of the central scatterer 
decreases. For the Sinai well a finite horizon, which means g(t)dt decreasing exponentially or 
faster, is compatible with the existence of periodic orbits of hyperbolic type (with no open 
corridor in the extended Lorentz version), whereas an infinite horizon, g(t)dt with algebraic 
tail, implies the existence of parabolic non-isolated periodic orbits that appear for R < R c . 
and determine the corresponding corridors. 

In the former case the decay law is exponential, and in the latter the decay law shows a 
crossover between an exponential and a power law decay (~ 1/t) for long times. 



IV. THE DECAY LAW IN D DIMENSIONS 

A. NUMERICAL STUDY OF THE DECAY 

In this subsection we show the results of the numerical study of the decay. It is a 
well known fact that, for short times, the behavior of the decay law in D dimensions is of 
exponential type ||. 

Our main purpose is to understand the behavior for long times. In other words, we study 
the long time tail of the decay law in D dimensions. As we will see in the following section, 



in order to extract information about the long time tail of the decay it is enough to study 
the function g(t)dt , the fraction of initial conditions for which the first collision with the 
scatterer occurs between t and t + dt, for the D dimensional system (to be more precise we 
will compute the integral of g(t), G = J t ,~ g(t') dt' ). From the numerical point of view, the 
study of G as a function of t instead of the decay law N(t)/N has a great advantage because 
it requires less CPU time and we can get better statistics with more initial conditions. We 
begin with Nq = 10 7 particles with random initial conditions in the phase space and the 
ratio Vq/E = 20. Fig. [j](a) shows the results of the numerical study of G for D = 2,3,4 
and radius of the D dimensional scatterer R = 0.23. The behavior is exponential for short 
times and becomes algebraic for long times (l/t s ). Fig. [l|(b) shows the tails of the G for 
D = 2,3,4 together with the best fit that predicts for all the curves an exponent 5=1. 
Fig. H|(a) is the same as Fig. |l](a) but for R = 0.4. For D = 2 the G is of exponential type. 
This agrees with the result of paper I. For the two-dimensional system R > R c \ = a/2/4 is 
compatible with a g(t) of "finite" horizon and the decay is exponential for all times. For 
D = 3, 4 again the best fit of the long time tail predicts an exponent 5=1. Fig. 0(b) shows 
this fit together with the numerical results of the long time tails for D = 3,4. These results 
suggest that for D > 2 there does not exist a critical value of the radius that changes the 
behavior of the decay for long times as happens for D = 2 (see I). In terms of the invariant 
trapped set, we can stress that for D > 2 there are always periodic orbits of parabolic type, 
and the initial conditions asymptotic to them are the ones that contribute to the long time 
tail of the decay. 

B. THEORETICAL STUDY OF THE DECAY 

One of the main results of the present paper is that for the Sinai well in D dimensions 
the long time tail of the decay law is (~ 1/t), independent of the dimension D. To our 
knowledge this is the first report in which this result is stated analytically and confirmed by 
numerical simulations. All the previous works only study the model in two dimensions [|7|], 



|| or conjecture for the D > 2 system an exponential decay of the velocity autocorrelation 
function for long times |J. To begin the theoretical study, we extend to the D dimensional 
case some results of the previous paper I that allow us to relate the decay law to the internal 
dynamics: 

Q(s) = U ° ° {S)/S ~ , (4.1) 

1 - (1 - UD )f( 8 ) 

where Q(s) = L[Q(t)] means the Laplace Transform and 

Q(t) = 1 - ^ (4.2) 

and N(t)/No is the fraction of particles present in the well at time t (the decay law). Taking 
into account the equation, 

dg/dt = -g(0)f(t) , (4.3) 

so that 

f(s) = 1 - sg(s)/g(0) , (4.4) 

we finally find 

A/ x _ up g(s)/s 

Wlj l + (l-u D )[g(s)s/g(0)-l]- ^ b) 

The preceding equation is a straightforward extension of the one that we have obtained in I 
for the two-dimensional problem, being in that case ujd=2 = w. To know Q(t), we must be 
able to inverse-transform ( f4.5|) . Computing the leading term of (|4.5| ) and employing ( |4.2|) 
we obtain 

N(t) ~u D [ g(t')dt' (4.6) 

Jt'=o 

We must center our attention on the function g(t)dt for our D dimensional problem to study 

the long time behavior ( [4.6[ ) of the decay law. 

We will extend the definition of corridors to higher dimensions. For that we must appeal 

to the Periodic Lorentz Gas in the .D-dimensional case which statistical properties in the 



hyperbolic domain (finite horizon) have been studied in detail in a recent reference || and 
a numerical study of some universal properties was performed in ||. We will say that a 
.D-dimensional Periodic Lorentz Gas has infinite horizon if the length of the free motion is 
unbounded. The first trivial extension is to define the corridors in D dimensions by the 
directions of the velocity v that satisfy Vi/vj = z\jz 2 V(z, j), where z\ and z 2 are coprimes 
integers. These are not the only directions that lead in D dimensions to unbounded free 
motion. For example for D = 3 the direction vi = 0, v 2 and V3 arbitraries (compatible 
with the condition that the modulus of the velocity, | v \, is one) is also a corridor in the 
meaning of free unbounded motion. The main characteristic to note is that as the dimension 
increases it is possible to find ever more directions that define unbounded motion for any 
value of the radius of the scatterer R < 0.5 which means that there is no critical value R c 
that destroys the corridors, as there was in two dimensions. However R c influences on the 
number of open corridors, since there are some corridors that disappear when R > R c . 
All the trajectories that contribute to the long time tail of the decay lie almost entirely in 
some corridor. In terms of the Sinai Well we say that these initial conditions are asymptotic 
to the parabolic periodic orbits that exist in D dimensions, leaving the bounded region with 
probability ujd after one collision with the D dimensional scatterer. 

At this point it is necessary to emphasize that apart from the corridors which asymptotic 
initial conditions contribute to the algebraic long time tail of the decay with an exponent 
5 = 1 (we call them principal corridors), there are other corridors in which initial conditions 
asymptotic to them lead to an algebraic tail of the decay (l/t M ) with the exponent \x > 1. 
A proof of the existence of these other corridors is given in Fig. 0. It shows for D = 4 and 
R = 0.4 two curves G as a function of t. In one (dashed line) all the corridors have been 
populated and the exponent for the long time tail is S = 1, in the other (solid line) we do not 
consider initial conditions asymptotic to the principal corridors, resulting in an exponent for 
the long time tail different from 1. 

Since the global behavior of the decay law for long times is (1/t) independent of the dimen- 
sion D, we call these other corridors hidden corridors, because there is no evidence of them 

8 



in the long time tail of the decay law. On the other hand the intermediate behavior of the 
decay law, that it is related to asymptotic conditions to the hidden corridors, results in a 
superposition of algebraic decays with exponents greater than and different from one, this 
behavior being more and more complicated as the dimension increases. 
We devote the remainder of this section to derive the explicit dependence on t of the dis- 
tribution g(t)dt for initial conditions that are asymptotic to the principal corridors using 
the D- dimensional Periodic Lorentz Gas model. Without loss of generality we will compute 
g(t)dt for the initial conditions that lie in the principal corridor (D = 3) defined by the 
directions of the velocity 1)1/1)2 = ±1 and 1*3 arbitrary that satisfy the condition | v \= 1, 
corresponding this case to R < R c . In fig. |] we show a schematic representation of the 
mentioned corridor of width I in which we have cut with the plane z = const, the spherical 
scatterers at their centers in order to simplify the figure. 

We remark that as | v \= 1 the distributions in times are equivalent to distributions in 
lengths. 

Let n(t*) be the fraction of initial conditions for which the first collision with some scatterer 
occurs in times t > t*. As we are interested in the large t behavior the angle a is proportional 
to n(t*) 

n(t*) ~ a ~ i (4.7) 

and 

- dn/dt* ~L^ g {t*) (4. 8 ) 

where g(t*)dt* is the fraction of initial conditions for which the first collision with the 
scatterer occurs between t* and t* + dt*. To be more precise we must consider the integration 
over the solid angle, but this gives 2tc, so we can conclude that inside a principal corridor 

git) - c f- (4.9) 



Use of the last expression in ( O ) leads to the mentioned behavior of the long time tail 
of the decay 



N(t) ~ — . (4.10) 

We stress that the last derivation is also valid for the other principal corridors. In the case 
of the hidden corridors, vestiges of some principal corridors when R > R c , the integration 
over the solid angle leads to an additional 1/t dependence, so in that case 

g(t) - |r ■ ( 4 - n ) 

The preceding arguments can be extended to higher dimensions. For the principal cor- 
ridors the dependence on t is exactly the same (~ 1/t). For the hidden corridors as the 
dimension D increases the integration over the solid angle gives rise to contributions of the 
type 1/t" with the 2 < fi < D. 

C. The Decay Law in terms of the Laplace Transform 

We begin computing the decay law N(t) as the Inverse Laplace Transform of a function 

g(A), 

rod 

N(t)= I q(\)e- xt d\ (4.12) 

where q(X) can be seen as the fraction of initial conditions that decay at time t and which 
decay rate varies between A and A + dX. 

In the context of ( |4.12| ) the whole decay law can be thought as the result of infinitely many 
decay processes of exponential 

type each one characterized by a decay rate A, this interpretation being independent of 
the dimension D. Our purpose is to relate the function q(X) with g(s) the Laplace Trans- 
form of the function g(t) which, as we have shown in the preceding subsection, determines 
univocally the decay law. Employing fl4.2|) and ( |4.12|) 



/oo 
q(X) e- xt dX (4.13) 

and making the Laplace Transform (LT) of ( [4.13D 
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^±d\= 1 --Q(s) (4.14) 

1° A + s s 

The left hand side of ( |4.14| ) is the Stieltjes Transform (ST) of the function g(A) and 



arises naturally from the iteration of the LT. Inverse transforming ( |4.14| ) we can determine 
univocally q(X) from Q(s). Here we do not give the exact expression for the inversion of 
fl4.14j ). For further details see [|IJ. We would like to emphasize that this point of view 
enables us to determine the fraction of initial conditions that decay at time t with a decay 
rate between A and A + dX when the function g(s) is known. 

V. SUMMARY AND CONCLUSIONS 

In the present work we have studied the decay of the Sinai well in D dimensions, this 
system being an extension to higher dimensions of the previous one studied in detail in I. The 
main difference between the D > 2 dimensional system and the analogue in two dimensions 
is that being both completely chaotic, the first one has an invariant parabolic trapped set for 
all the values of the radii R of the scatterer, whereas in the second one the invariant trapped 
set can be fully hyperbolic or have a parabolic trapped subset if R < ^. In terms of the 
decay law it means that for the present system the decay law is always of algebraic type for 
long times (related to initial conditions asymptotic to the parabolic periodic orbits in D > 2 
dimensions). We have related the decay law in D dimensions to internal distributions that 
characterize the dynamics and concluded that the exponent 5 for the algebraic long time 
tail is S = 1 irrespective of the dimension. This behavior is in agreement with the results 



encountered in ref . JTTJ • In that work the authors investigate in D dimensions a system like 
a Sinai billiard and allow the decay providing a small window in one of the walls of the 
container. They conclude that for the regular case (without a scatterer in the center) the 
decay law is algebraic (~ 1/t) independent of D. As the trapped set of the regular system 
coincides with the parabolic periodic orbits of the Sinai well, the result is consistent with 
our conclusion. 
We also studied the intermediate behavior of the decay law and concluded that it is related 

11 



to the so called hidden corridors, which asymptotic initial conditions contribute to the decay 

law with a temporal dependence of the type (l/i M ), with 2 < yU < D. 

The present work was partially supported by CONICET and Fundacion Antorchas. 
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APPENDIX A: 

In this appendix we derive the explicit expression for ud, the transition probability from 
the bounded to the free region after one collision with the scatterer. It can be evaluated 
using ergodic theory as the ratio between all orientations of momentum in the free region 
and all possible orientations of momentum. 

The space of momenta is a D dimensional unit sphere (the modulus of p is one). Let 
d^um{D) be the solid angle subtended by \I/« m , where ^u m is given in expression (|2.2|) , and 
Q(D) the total solid angle in D dimensions. As we have established in (|2.3| ) 



UD = ^m~ ' (A1) 

so for computing up we must have the explicit expressions for Q(D) and d^>u m (D). 
Introducing the polar coordinates in D dimensions (D > 3) 

(r,(f>,8i, ,0d- 2 ) (A2) 

with r = 1, < < 2tt and < 9i < it, i = 1, ....D — 2 we can derive the expression 

D-2 

Q(D) = 2tt [] / sin 1 ^ d9 { . (A3) 

Employing the well-known formula [12| 

f sin 4 e t d6 t = v^F r j' + l,2 \ (A4) 

Jo V T(i + 2/2) v ; 

and putting it back in ( |A3| ) , we obtain 

n D - 2 
n(D) = — — -- (A5) 

v ; rp/2) v ; 

In order to compute d^fu m (D) we must perform the integration in the variables 9i with 
i — 1, ...D — 3 between < 9i < 7r, and < 9d-i < ^iim 

<M lim (D) = 2tt J] / sm^idOi sin D - 2 9 D _ 2 d9 D _ 2 (A6) 

i=1 Jo Jo 

after a straightfoward calculation we obtain 
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and from (|2.3| ), 



d^um(D) 



2 (v^r- 1 * 



0-1 iTr-D-l 

Urn 



and 



Wfl 



(£> - 1) r(D - 1/2) 
results 

2 rp/2) Eg- 1 

0F (£> - 1) T(D - 1/2) ' 



(A7) 



(A8) 
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FIGURES 
FIG. 1. (a)Numerical results of G for D = 2, 3, 4. We have fixed a (the side of the D dimensional 



square well) as the unity of length and a/\/2(E + Vq) as the unity of time. The log-log plot shows 
G vs. t for radius R = 0.23 of the central scatterer. (b)The long time tail together with the best 
fit to G vs. t that is consistent with the exponent 5 = 1 for D = 2,3, 4. 

FIG. 2. (a)Numerical results of G for D = 2,3, 4. The units are the same as in FIG.|l[ The 
log-log plot shows G vs. t for radius R = 0.4 of the central scatterer. (b)The long time tail together 
with the best fit to G vs. t that is consistent with the exponent 5 = 1 for D = 3,4. 

FIG. 3. Numerical results of G for D = 4. The units are the same as in FIGJI|. The log-log 
plot shows G vs. t for radius R = 0.4 of the central scatterer. The dashed curve corresponds 
to initial conditions in all phase space. The solid curve results from having not populated any 
principal corridors. The algebraic long time tail is originated by the hidden corridors. The best fit 
is consistent with an exponent 5 = 2. 

FIG. 4. Periodic configuration of scatterers in D = 3 dimensions, for a constant value of the 
coordinate z. The coordinate z has been chosen in order to cut the spherical scatterers at their 
centers. The principal corridor defined by v\/v2 = 1 and V3 arbitrary that satisfied the condition 
I if |= 1 is shown. The width of the corridor is /, and the angle a (~ l/t*) is proportional to the 
number of initial conditions that collides with some scatterer in a time t > t* . 
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